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ABSTRACT 


MSTRIP  has  undergone  a  succession  of  improvements 
since  its  original  publication  in  1968.  Perfor¬ 
mance  of  the  1974  and  later  versions  relating  to 
aspects  including  weakly  coupled  (<-30  dB)pairs, 
low  shield  heights  (closer  than  substrate  thick¬ 
ness),  and  extremely  low-impedance  lines  will  be 
discussed.  A  FORTRAN  list  of  the  1978  version  is 
appended . 
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MSTRIP2 :  PARAMETERS  OF  MICROSTRIP  TRANSMISSION 
LINES  AND  OF  COUPLED  PAIRS  OF  LINES  - 
1978  VERSION  AND  ITS  APPLICATION 

I.  INTRODUCTION 

The  "quasi- numerical  analysis  of  the  normal  modes  of  propagation  on 
coupled  pairs  of  luicrostrip  transmission  lines  named  MSTRIP*  published1  by 
T.G.  Bryant  and  J.A.  Weiss  in  1968,  incorporated  a  Green' s-function  represen¬ 
tation  of  "bound  charge"  at  the  upper  surface  of  the  dielectric  substrate. 

For  a  wide  range  of  applications,  generally  those  in  which  the  cross- 
sectional  dimensions  are  small  enough  relative  to  wavelength,  the  accuracy  of 
this  method  is  very  good;  in  fact,  results  of  MSTRIP  have  been  widely  cited 
over  the  years  as  a  standard  for  assessment  of  the  accuracy  of  other  micro- 
strip  algorithms.  MSTRIP  has  been  used,  apparently  without  modification,  in 
commercially  available  linear  circuit  analysis  packages.2 

The  principal  approximations  (other  than  the  usual  substitution  of  dis¬ 
crete  in  place  of  certain  continuous  variables)  are  (a)  the  quasi-TEM  assump¬ 
tion,  whereby  the  results  of  the  analysis  are  all  based  on  the  determination 
of  static  capacitance,  and  (b)  the  assumption  of  zero  thickness  and  perfect 
conductivity  of  the  conductors.  Numerous  papers  have  appeared  in  which  this 
basic  analysis  has  been  adopted  as  the  starting  point  for  determination  of 
the  influences  of  dispersive  effects  and  of  finite  conductivity  and  thickness 
of  the  strips. 

The  program  has  had  some  revisions  for  improved  speed,  accuracy,  and 
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range  of  applicability  during  the  years  since  its  original  publication.  We 
present  this  report  in  order  to  bring  to  the  attention  of  microstrip  circuit 
users  the  current  version  of  MSTRIP  and  its  capabilities  and  to  address  some 
questions  which  have  arisen  in  relation  to  recent  application  requirements. 

In  1970  and  1971,  an  improved  program  was  formulated3*4  in  which  the 
tabulated  dielectric  Green's  function  was  replaced  by  a  Fourier-integral 
evaluation.  This  change  resulted  in  a  very  substantial  improvement  in  speed 
and  simplicity,  better  accuracy,  and  removal  of  the  limitations  on  the  range 
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of  substrate  dielectric  constant  values.  It  also  introduced  the  option  of 
including  the  influence  of  a  shield,  or  upper  ground  plane,  above  the  circuit 
surface.  The  geometry  of  the  shielded,  coupled  pair  of  microstrip  lines  is 
illustrated  in  Fig.l.  A  detailed  review  of  methods  and  results,  with  a 
FORTRAN  list  of  the  version  of  the  program  then  current,  was  published5  in 
1974.  Following  some  improvements  made  in  1978,  the  program  name  was 
modified  to  MSTRIP2.  The  FORTRAN  source  of  this  program  is  shown  in  the 
Appendix. 

In  more  recent  times,  users  have  brought  up  questions  relating  to  the 
applicability  and  accuracy  of  MSTRTP  for  specific  application  requirements. 
Examples:  the  questions  of  upper  limits  on  substrate  dielectric  constant  and 
on  strip  widths;  usability  of  the  program  for  shielded  aicrostrip  with  the 


Fig.l.  Shielded  microstrip  structure. 
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shield  very  close  cc  the  substrate  surface  (H2/H1  <‘2);  the  accuracy  of  the 
program  in  cases  of  very  weakly  coupled  pairs  of  lines,  for  which  the  differ¬ 
ences  between  the  two  normal-mode  characteristic  impedances  and  velocities, 
which  determine  the  backward-  and  forward-wave  coupling  coefficients,  are 
very  small. 

Applications  of  these  unusual  microstrip  structures  have  arisen  in  vari¬ 
ous  contexts.  For  example,  low-impedance  (hence  wide-line)  transmission 
lines  are  needed  to  match  the  low  impedance  of  Josephson- junction  logic 
elements  in  superconducting  computers.6  Low  shield  heights  can  be  used  to 
correct  the  even-  and  odd-mode  velocity  differences  of  coupled  lines  on 
anisotropic  dielectrics  such  as  sapphire.7  Weakly  coupled  pairs  of  lines  are 
used  as  directional  couplers  in  large  time-bandwidth  dispersive  filters.8 

In  the  following  sections,  information  is  presented  relating  to  such 
structures  and  to  some  tests  of  accuracy  of  the  program  which  their  consider¬ 
ation  provides. 

II.  SUBSTRATS  DIELECTRIC  CONSTANT 

Nothing  in  MSTRIP2  imposes  a  limit  on  the  substrate  dielectric  constant 
K.  The  program  has  served  successfully  with  K  values  from  1  to  6000  (in 
fact,  values  of  K  less  than  unity  have  been  used  in  order  to  represent  an 
inverted  structure  -  see  Sec. HI). 

III.  APPLICATION  TO  CASES  OF  A  CLOSELY  SPACED  SHIELD 

For  values  of  %/%  less  than  2  (see  Fig.l),  the  gap  between  the  shield 

and  circuit  surface  is  smaller  than  the  substrate  thickness.  In  such  cases, 

MSTRIP2  loses  accuracy.  For  example,  with  *  1.25,  S/H^  *  3.5,  W/Hj  * 

0.50,  and  K  *  9.6  (alumina),  MSTRIP2  gives  the  physically  unreasonable  result 

that  Z  <  Z  . 
oe  00 
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TABLE  I 


TRANSFORMATION  FROM  STANDARD  TO  INVERTED  COMPUTATION 
IN  THE  CASE  OF  A  CLOSELY  SPACED  SHIELD* 


Unprimed:  Parameters  of  interest  as  used  in  the  standard 

computation  (Fig.l). 


Primed: 


Input/output  parameters  to  be  used  to  analyze  the 
some  structure  by  inverted  computation. 


MSTRIP  Input:  Hj  -  H2  -  Hj 


/wv.  VJ1 1_ 

\H j/  R  -  1 


R  -  1 


[nj  r  -  1 


r’«  1 


To  interpret  MSTRIP  output,  transform  tabulated  values  as  follows. 
(Primed  variables  are  output  of  the  inverted  computation,  unprimed 
variables  are  those  of  the  actual  structure  of  interest.) 


Keff  *  Kiff* 


C  -  C’K 


*See  Sec. III. 


This  deficiency  may  be  circumvented  b'’  inverting  the  problem,  placing 
the  gap  (K  *  1)  on  the  bottom,  and  the  substrate  (K  >  1)  above,  so  that 
H2/Hi  >  2.  Because  MSTRIP  incorporates  the  assumption  that  the  gap  has  a 
dielectric  constant  of  unity,  the  input  and  output  parameters  must  be  trans¬ 
formed  according  to  the  system  shown  in  Table  I. 

Using  this  procedure  as  a  check,  it  may  be  verified  that  MSTRIP2  is 
accurate  to  values  less  than  1.5.  For  example,  for  a  50-fl  pair  of 

coupled  lines  on  sapphire  with  =  1.5,  S/Hj  *  3.5,  the  two  methods  agree 

to  within  0.01  percent  in  impedances  and  0.02  percent  in  velocities.  The  two 
agree  to  within  0.04  dB  in  backward  coupling,  even  though  it  is  a  weak  -64  dB. 
With  the  shielding  height  reduced  to  H2/H1  *  1.25,  the  discrepancy  in  imped¬ 
ance  between  the  standard  computation  and  the  inverted  structure  increases  to 
0.3  percent,  unacceptable  for  the  calculation  of  coupling  in  this  weakly 
coupled  structure. 

IV.  COMPARISON  WITH  ANALYTIC  CALCULATION 

In  the  special  case  H2/H^  *  2,  the  quasi-static  electric  field  in  the 
shielded  (coupled)  microstrip  structure  with  zero  conductor  thickness  may  be 
shown  to  be  symmetric  about  the  dielectric/air  interface.  Consequently,  the 
quasi-static  characteristics  of  the  microstrip  (MS)  are  identical  to  those  of 
a  balanced  stripline  (SL)  structure  of  the  same  dimensions  filled  with  a 
material  of  dielectric  constant  Kg^  »  (1/2)  (K^g  +  1). 

This  equivalency  provides  a  convenient  check  on  shielded  microstrip 
calculations,  as  analytical  expressions  for  coupled  striplines  are  well 
known.  In  Table  II  are  shown  the  results  of  calculations  for  several  values 
of  S/H^  and  W/Hj,  using  MSTRIP  and  MSTRIP2  for  shielded  microstrip  (K^*  9.6) 
and  an  analytic  approximation  for  the  stripline  (Kg^  *  5.3).  A  transcenden¬ 
tal  approximation  valid  to  8  parts  in  10^  compared  to  the  exact  elliptic 
integral  expression  for  stripline  impedances  was  used,  given  as  Eqs.6.3.1, 
6.3.2,  3.2.7,  and  3.2.8  in  the  book  by  Gunstoh.9  Both  MS"RIP  and  MSTRIP2 
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TABLE  II 


IMPEDANCE,  COUPL] 
CALCULATED  BY  MSTRIP,  1 
K  -  9.1 

[NO,  AND  EFFECTIVE  DIELECTRK 
1STRIP2 ,  AND  BALANCED  STRIPL1 
FOR  MiCROSTRIP,  5.3  FOR  STf 

CONSTANT  VALUES 

NE  METHODS,  H2/H!  - 
tIPLINE 

2; 

Z  (Q) 

Z  (0) 

-k  (dB) 

oe 

oo 

B 

Strip- 

Strip- 

Strip- 

(K  cJ 
eff  even 

S/H! 

W/Hj 

MSTRIP 

MSTRIP2 

line 

MSTRIP 

HSTRIP2 

line 

MSTRIP 

MSTRIP2 

line 

(MSTRIP) 

0.5 

0.3 

90.29 

89.82 

89.36 

58.26 

53.26 

57.77 

13.3 

13.4 

13.4 

5.346 

0.8 

57.72 

57.25 

56.84 

*0.60 

40.60 

40.07 

15.2 

15.4 

15.2 

5.372 

1.3 

42.94 

42.46 

42.12 

32.43 

32.43 

31.90 

17.1 

17.5 

17.2 

5.397 

1.0 

0.3 

81.75 

81.27 

80.81 

67.29 

67.29 

66.81 

20.3 

20.5 

20.5 

5.351 

0.8 

53.55 

53.07 

52.65 

45.63 

45.63 

45.16 

22.0 

22.5 

22.3 

5.378 

1.3 

40.55 

40.07 

39.70 

35.60 

35.60 

35.16 

23.7 

24.6 

24.3 

5.403 

2.0 

0.3 

76.25 

75.77 

75.31 

72.88 

72.89 

72.42 

32.9 

34.2 

34.2 

5.354 

0.8 

50.69 

50.22 

49.78 

48.68 

48.68 

48.23 

,  33.9 

36.2 

36.0 

5.382 

1.3 

38.87 

38.39 

38.00 

37.47 

37.47 

37.06 

34.7. 

38.3 

38.1 

5.407 

4.0 

0.3 

74.87 

74.39 

73.93 

74.27 

74.27 

73.80 

47.9 

61 .5 

61.5 

5.355 

0.8 

49.96 

49.48 

49.05 

49.42 

49.42 

48.98 

45.2 

63.5 

63.3 

5. 383 

1.3 

38.43 

37.95 

37.55 

37.91 

37.91 

37.51 

43.4 

65.6 

65.3 

5.408 

10. 0 

0.3 

74.81 

74.33 

73.87 

74.33 

74.33 

73.87 

49.9 

>100 

>100 

5.355 

0.8 

49.93 

49.45 

49.01 

49.45 

49.45 

49.01 

46.4 

>100 

>100 

5.383 

1.3 

38.41 

37.93 

37.53 

37.93 

37.93 

37.53 

44.1 

>100 

>100 

5.408 

were  run  with  M  *  20  substrips  per  strip.  In  addition  to  the  even-  and  odd¬ 
mode  impedances,  the  backward-wave  coupling  strength  for  a  quarter-wave 
length  of  line  is  shown,  given  by  the  expression 


h  (dB)  “  20  108io 


oe 


oo 


The  even-mode  effective  dielectric  constant  given  by  MSTRIP  is  also  shown. 

The  odd-mode  constant  given  by  MSTRIP  and  both  mode  constants  given  by 
MSTRIP2  exhibit  the  proper  value  of  5.3. 

It  is  evident  that  both  versions  of  MSTRIP  give  reliable  impedance  and 
backward  coupling  values  for  k^  in  the  range  of  stronger  coupling  (S/H^  <  1, 
k_  >  -25  dB).  Impedances  are  within  1  8  of  the  stripline  values  and  coupling 

O 

strengths  are  within  1  dB.  In  the  -30-  to  -40-dB  range  (S/Hj  ~  2)  MSTRIP2 
retains  its  accuracy  in  coupling  strength,  while  errors  of  more  than  3  dB 
appear  in  MSTRIP.  For  S/Hj^  >  2,  MSTRIP  rapidly  loses  accuracy,  never  pre¬ 
dicting  coupling  weaker  than  -50  dB.  MSTRIP2  remains  within  a  few  tenths  of 
one  decibel  of  the  stripline  value. 

It  is  noteworthy  that,  over  the  whole  range  of  S/H^  and  W/Hj  tested,  the 
odd-mode  impedances  given  by  MSTRIP  and  MSTRIP2  are  identical  and  exceed  the 
strip line  values  by  0.4  to  0.5  fl.  MSTRIP2  gives  ZQe  values  which  are  simi¬ 
larly  0.4  to  0.5  0  higher  than  stripline  values,  while  MSTRIP  gives  Z 

oe 

values  about  1  0  too  high.  Thus,  even-  and  odd-mode  impedance  errors  given 
by  MSTRIP2  track  each  other,  giving  accurate  coupling  values,  while  a  roughly 
0.5-0  difference  exists  between  the  MSTRIP  even-  and  odd-mode  impedance 
errors,  resulting  in  large  errors  in  coupling  strength  for  weakly  coupled 
lines. 

There  is  also  a  significant  error  in  the  even-mode  effective  dielectric 
constant  given  by  MSTRIP.  This  error  is  positive,  opposite  in  sign  to  that 
expected  on  the  basis  of  the  ZQg  error.  The  error  results  in  a  predicted 
even-mode  phase  velocity  less  than  that  of  the  odd  mode.  On  this  basis,  a 
forward  coupling  strength  of  -40  to  -50  dB  over  a  90*  line  segment  would  be 
predicted  for  all  lines  considered  in  Table  I. 
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MSTRIP  is  adequate  for  design  of  couplers  tighter  than  -25  dB,  although 
it  may  give  inaccurate  values  of  directivity.  For  calculation  of  cross  talk 
between  closely  spaced  lines  in  applications  such  as  MICs  where  space  is  at  a 
premium,  the  current,  more  accurate  version,  MSTRIP2,  is  needed. 

V.  APPLICATION  TO  EXTREMELY  LOW- IMPEDANCE  LINES 

In  1980  an  MSTRIP  user,  J.C.  Cozzie,  pointed  out1®  a  flaw  which  he  had 
discovered  when  he  attempted  to  determine  microstrip  parameters  for  extremely 
low-impedance  lines:  Zq  <  3  Q,  which  requires  strip  widths  of  W/Hj^  >  70  in 
the  case  K  *  2.45  which  was  of  interest  to  him.  Cozzie  located  the  flaw  in 
the  subroutine  MGREEN,  recognized  it  as  a  consequence  of  an  approximation 
embedded  in  the  integration  procedure,  and  described  how  he  eliminated  it 
with  changes  including  the  substitution  of  a  Gaussian  quadrature  procedure. 
(Cozzie  was  working  with  the  1971  or  1974  version  of  MSTRIP.  By  coincidence 
we  had  also  converted  MGREEN  to  Gaussian  quadrature  integration,  calling  IBM 
Scientific  Subroutine  DQG32,  in  the  1978  version.) 

At  Cozzie* s  instigation,  we  modified  certain  details  of  MGREEN  to  remove 
the  flaw.  At  the  time,  we  had  the  following  further  thoughts  regarding  very 
wide  strips.  First:  if  the  scale  of  dimensions  is  such  that  the  width  be¬ 
comes  comparable  to  a  half-wave length  (in  a  medium  of  dielectric  constant 
Keff  as  evaluated  by  MSTRIP2)  in  the  frequency  range  of  interest,  then 
spurious  propagation  effects  and  coupling  problems  are  to  be  expected.  In 
that  case,  it  would  probably  be  advisable  to  seek  some  alternative  means  to 

obtain  the  circuit  function  contemplated.  Second:  for  widths  of  W/Hi  >  20 
or  so.  very  accurate  results  can  be  obtained  by  decomposing  the  strip  capaci¬ 
tance  into  the  sum  of  W/Hi  times  a  parallel-plate  part  per  unit  increment  of 
W/Hj  plus  an  edge  contribution.  Once  suitable  limiting  values  of  those  two 
capacitance  parameters  have  been  determined  by  a  few  trials,  the  microstrip 
parameters  for  all  cases  of  wider  strips  with  the  same  substrate  K  value  can 
be  quickly  determined.  Cozzie  noted  this  suggestion  and  briefly  presented 
his  version  of  it  in  his  paper.10 
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-M5TRIP2- 

PARAMBTBRS  OP  MICROSTRIP  TRANSMISSION  LIHES  AMD 
OP  COUPLED  PAIRS  OP  MICROSTRIP  LINES 

IMPROVED  —  VERSION  OP  ADO.  1  ,  1978 

RBPSR8NCSS  -  1.  T.Q.  BRYANT  AND  J.A.  WEISS,  IEEE  TRANS  MTT-16  1021 

(DECEMBER  1968). 

2.  J.A.  WEISS  AND  T.Q.  BRYANT,  BLBCTR03ICS  LETTERS 

VOL.  5  P.  517  (OCTOBER  IS,  1S69). 

3.  J.A.  WEISS  AND  T.O.  BRYANT,  ELECTRONICS  LETTERS 

VOL.  6  P.  462  (JULY  23,  1970) ;  SEE  ALSO  ERRATUM, 

ELECTRONICS  LETTERS  VOL.  6  P.  5S0  (AOS0ST  20,1970). 

4.  T.3.  ERYANT  AND  J.A.  WEISS,  IEEE  TRANS  MTT-19  418 

(APRIL,  1971). 

5.  J.A.  WEISS,  ADVANCES  IN  MICRtWAVBS,  VOL.  8, 

PP.  295-320;  ACADEMIC  PRBSS,  1974. 

FOR  MICROSTRIP  COMPRISING  -  A  DIELECTRIC  SUBSTRATE  OP  RELATIVE 
DIELECTRIC  PERMITTIVITY  X  AND  THICKNESS  HI  LYING  ON  A  CONDUCTING 
GROUND  PLANE;  A  SINGLE  STRIP  OF  WIDTH  W  OR  TWO  PARALLEL  STRIPS  OP 
EQUAL  WIDTH  W  WITH  INNER  EDGES  DISTANCE  S  APART,  LYING  OK  THE  UPPER 
SURFACE  OF  THE  SUBSTRATE;  THICKNESS  OF  THE  STRIP  MATERIAL  ASSUME) 

TO  BE  NEGLIGIBLE;  AH  UPPER  GROUND  PLANE  LOCATED  PARALLEL  TO  AND 
DISTAHOB  H2  ABOVE  THE  LOWER  GROUND  PLANE  (OPTIONAL) . 

THE  PROGRAM  COMPUTES  THE  FOLLOWING  TRANSMISSION-LINE  PARAMETERS 
FOR  A  SINGLE  STRIP  OR  FOR  BOTH  THE  EVEN  AND  ODD  MODES  OF  A  COUPLED 
PAIP  OF  STRIPS  - 

CHARACTERISTIC  IMPEDANCE; 

PHASE  VELOCITY; 

EFFECTIVE  DIELECTRIC  CONSTANT  K-SFFf 

CAPACITANCE  PER  UNIT  LENGTH  PER  STRIP  (OPTIONAL). 

THE  CALCULATION  BffiODIBS  THE  QUASI-STATIC  APPROXIMATION,  AS 
EXPLAINED  IN  RBFBRENCB  1;  THE  APPROXIMATION  IS  ACCURATE  PROVIDED 
HI  IS  NO  GREATER  THAN  A  FEW  PERCENT  OF  THB  WAVELENGTH  IB  JUi  UN¬ 
BOUNDED  MEDIUM  OP  PERMITTIVITY  K. 

THB  PROGRAM  ACCBPTS  ANY  VALUB  OF  X,  S/HI,  AND  W/HJ . 

INPUT  AND  OUTPUT  MEDIA  ARB  IDENTIFIED  IN  THB  PROGRAM  BY  THE 
FOLLOWING  NAMES  -  INPUT  (USUALLY  A  CARD  READER),  IV;  OUTPUT  (USU¬ 
ALLY  A  LINS  PP.1NTBR) ,  IW. 


•MST00010 

MST00020 

MST00030 

MST00040 

MST00050 

MST00060 

MST00070 

MST00080 

MST00090 

MST00100 

M6T00U0 

MST00120 

MSTC'0130 

MST00140 

H3T00150 

MBT00160 

M6TS0170 

MSTG0180 

MST00190 

MST00200 

MST00210 

KST00220 

MST00230 

MST00240 

M6T0025P 

MST00260 

MST00270 

MST00280 

MGT00290 

MST00300 

MST00310 

M6T00320 

MST00330 

MST00340 

KST00350 

MST00360 

MST00370 

M6T00380 

HST00U9Q 

MSTJIHOC- 

M5T00410 

MGT0O4I0 

MOT00430 

MST00440 

MST00450 

MST00460 

MST00470 

MST00480 

MST09490 
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C  INPUT  DATA  -  THE  INPUT  DATA  CARD  IS  WRITTEN  IN  NAMELIST  FORMAT  MSTOObOO 
C  AND  CONTAINS  SIX  DECIMAL  NUMBERS  AND  ONE  INTEGER  NUMBER,  ALL  MST00510 
C  SEPARATED  BY  COMMAS,  AS  IN  THE  FOLLOWING  EXAMPLE.  MST00520 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


MST00530 

SCONST  WH1«0.1,DELW-').2,NT»20,R-0.0,DIEK»S.6,SH1»0.4,.AtR,i.o,  SEND  MSTO0S40 

MST00550 

THIS  CARD  SPECIFIES  THB  VALDES  OF  PROGRAM  VARIABLES  HAVING  THE  FOL-  MST00560 
LOWING  NAMES  -  MST00570 

HST00580 

WH1 . DELW,NT , R , DIKK ,SHI , AIR  MST00590 

MST00600 

WH1  IS  THE  STARTING  VALUE  OF  W/Hlj  DBLW  IS  THB  IHCRBMENT  OF  W/Hl;  MST00610 

NT  IS  THB  NUMBER  OF  IINBS  IN  THB  OUTPUT  TABLB  (EQUAL  TO  THB  HUM-  MST00620 

BER  Ob'  VALOBS  OF  W/Hl)  j  R  IS  THB  GROUND  PLAHB  RATIO  H2/H1  (FOR  THB  MST00630 

OPTION  OF  NO  UPPER  GROUND  PLAHB,  SET  R  -  0.0))  DISK  IS  THB  RELA-  MST00640 

TIVB  DIBLBCTRXC  PERMITTIVITY  K  OF  THB  SUBSTRATE  I  SHI  IS  THB  VAI.UB  MST00650 

OF  S/Hlj  AIR  TAXBS  ONE  OF  THB  TWO  VALUES  0.0  (SIBGLB  STRIP)  OR  MST00660 

1.0  (COUPLED  STRIPS).  THUS  IN  THB  EXAMPLE  THB  OUTPUT  TABLB  WILL  RS-  KST00670 

FBR  TO  THB  EVEN  AND  ODD  MODES  OF  COUPLED  STRIFE  WITH  SPACING  S/HI  MST00680 

-  0.4  ON  SUBSTRATE  OF  PERMITTIVITY  K  »  9.6,  HO  UPPER  GROUND  PLANE,  MST0D690 

AND  WILL  LIST  20  LIMBS,  FROM  W/Hl  >0.1  THROUGH  W/Hl  -  2.0  IN  MST00700 

STEPS  OP  0.2.  MST00710 

MST00720 

OUTPUT  -  THB  OUTPUT  FORMAT  IS  ILLUSTRATED  IB  THB  FOLLOWING  MST00730 

TABLB  SHOWING  THB  FIRST  LINE  OF  OUTPUT  DATA  FOR  THB  ABOVE  EXAMPLE.  MST00740 

MST00750 

HST00760 

H2/H1  -  0.0  K  -  9.60  COUPLED  STRIPS  6/HI  -  0.4  MST00770 

M6T00780 

M6T00790 

W/Hl  20(B)  ZO(O)  V (B)  V(O)  K-EFF(B)  K-B7F(0)  MST00800 

OHMS  H*08  M/SBC  M6T00810 

MST00820 

0.100  141.210  77.250  1.214  1.299  6.102  5.328  MST00830 

MST00840 

MST00850 

SUBROUTINES  REQUIRED  -  THE  PROGRAM  CALLS  THREE  SUBROUTINES  OB-  MST00860 
TAINBD  FROM  THE  IBM  SCIENTIFIC  SUBROUTINE  PACKAGE  (SSP).  THEY  ARE  MST00870 
DGBLS,  SICI,  AND  DQG32.  FOR  DESCRIPTIONS  AND  LISTS  OF  THESE  PROGRAMS  MST00880 
SBB  THB  IBM  SSP  MANUAL.  IN  ADDITION,  THB  PROGRAM  INCLUDES  FIVB  SUB-  MST00890 

ROUTINES  OP  ITS  OWN  -  OUTPUT,  XQBH,  AMAT,  MPHI,  AND  M9RBEN  AND  ONB  MST00900 

FUNCTION  SUBPROGRAM  -  GINT.  NOTE  ON  THB  IMPROVEMENT  DATED  AUG.  1,  MST00910 

1?7«,  THB  INTEGRATION  SUBROUTINE  DQG32  (32-POINT  GAUSSIAN  QUAD-  MST00920 

RATURB  METHOD)  REPLACES  7H3  SUBROUTINE  DQSF (SIMPSON* S  RULE)  ,  MST00930 

RESULTING  IN  ACCURACY  SUBSTANTIALLY  BETTER  THAN  CSS  PERCENT  IN  ALL  MST00940 
CASES  TESTED.  REVISIONS  AND  TBSTS  PERFORMED  BY  R.  C.  LEWIS.  MST0G550 

MST00960 

NOTE  ON  DIMENSION  REQUIREMENTS  -  FOR  TH3  DETERMINATION  OF  MST00970 

CHARGE  DISTRIBUTION  BACH  STRIP  IS  DIVIDED  INTO  M  SUBSTRIPS.  IK  MST00980 

THB  PRBSBNX  PROGRAM  M  »  20,  BUT  THB  PARAMBTBR  M  CAN  BB  INCREASED  KST00990 
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C  FOR  IMPROVED  ACCURACY  OR  REDUCED  FOR  IMPROVED  SPEED  OF  BJI2C0-  MSTOiOOC 
C  TION,  THE  DIMENSIONS  OF  THB  VARIABLES  DEPEND  ON  THB  VALCB  OP  M  MST01010 
C  AS  FOLLOWS  -  MST01020 
C  MST01030 
C  V (M) , AUX(M-l) ,X(3M-1) ,PHI(3M-1) ,A(M(M*1) /2) , B (M(M< 1) /2)  MST01040 


C 

C 

IMPLICIT  RRAL«8(A-H,0-Z) 

DIMENSION  V(20) ,AOX(19) ,GO(2) 

COfWON/DIM/X(59)  , PHI (59)  ,A(210j  ,B(210) 

NAM3LI ST/C0NST/WH1 , DRLW ,  NT ,  R ,  DI EK ,  SHI , AIR 
DATA  OO/ '  Y* , 'N'/ 

CALL  TIMBS  (DATB,TIMB) 

IV-5 

IW-1 

M-20 

H-l 

N?.iTB(IW,104) 

801  00  901  NRBP-1,25 
BPS-1. 08-07 

IF (NRBP  .BQ.  1)  OO  TO  91 
92  WRITB(6 ,102) 

READ(IV,103)  00  YN 

IF (QOYN  .BQ.  00(1))  00  TO  91 

IF (GOYN  . NB.  00(1)  .AND.  OOYK  .HE.  00(2))  00  TO  92 
GO  TO  999 
91  WRITS (6, 101) 

RBAD( IV, CONST) 

IF (AIR)  38,7,8 

7  WRITS (IW,  105)  R,DIBK 
WRITS (IW, 106) 

00  TO  802 

8  WRITS (IW,  107)  R, DISK, SHI 
WRITS (IW, 108) 

802  DO  902  K«1,ST 
AK-K 

WH-WH1+ (AX~ 1.0) «DBLW 
CALL  XQSN(M,NH,SH1) 

NH-AIR+1 

ADIKK-DIEK 

803  DO  903  IJ-1,2 
IF(IJ.EQ.l)  DIBK-1.0 
IF (IJ . BQ. 2)  DISK-ADIKK 

IP(R. BQ. 0.0. AND. DISK. BQ. 1.0)  00  TO  39 
CALL  M3RBKN  (M,WH,  SHI , DISK, R) 

00  TO  804 

39  CALL  MPHI (WH,M,ATR) 

804  DO  903  JJ»1,NN 
CALL  AMAKA1P..M) 

DO  23  I«1,M 


MST01050 

•***MST01060 

MST01070 

M6T01080 

MST01090 

MST01100 

MST01110 

MST01120 

MST01130 

MST01140 

MST01150 

MST01160 

HST0U7C 

M6T01180 

MST01190 

MST01200 

MST01210 

MST01220 

MST01230 

MST01240 

HST01250 

MST01260 

MGT01270 

MKE01280 

MST01290 

MfcTulICC 

M6T01310 

MST01320 

MST01330 

M6T01340 

M6T01350 

MST01360 

MST01370 

MST01380 

MST01390 

MST01400 

KST01410 

KST01420 

MST01430 

MST01440 

MST01450 

MST01460 

MST01470 

MST01480 

MST01490 

MST01500 
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23  V(I)=1.0  .  •  MST01510 

IF ( JJ  .BO.  1)  CALL  DGSLS{V,A,M,N,BPS,IER,AUX)  MST01520 

IF ( JJ  .EO.  2)  CALL  DGSLS ( V , B ,M,N , EPS , IER, AUX)  MST01530 

IF'IER.NB.O)  WRITE ( IW, 111)  IER, EPS  MST01540 

111  FORMAT  ( '  IER«  ',13,',  IN  SUBROUTINE  D3BLS,  SO  THE  CHARGB  DENSITY  MTT01550 
2  COULD  NOT  BB  CALCULATED  TO  THE  PRECISION  OP  ',813.6,*  DIGITS')  MST0156"! 


CAPSUM-0 . 0 
BO  21  I-1,M 

21  CAPSUM-CAPSUM«-V(I) 

CC-CAPSUM*111.256 

IF (JJ  .BO.  1  .AND.  IJ  .EO.  1)  CAP1E-CC 

IF ( JJ  .BQ.  1  .AND.  IJ  .BO.  2)  CAPXB-CC 

IPCJJ  ,SQ.  2  .AND.  IJ  .BQ.  1)  CAPlO-CC 

IF (JJ  .BO.  2  .AND.  IJ  .BQ.  2)  CAPXO-CC 

903  CONTINUE 

CALL  OUTPUT (CAP IB , CAPKB , CAPIO ,CAPKO ,WH , IW,NE) 

902  CONTINUB 

WRITS ( IW, 110 )  DATS , TIME 
WRITS (IW, 19 9) 

GO  TO  901 
38  WRITB (6 , 109) 

101  FORMAT ( / , 3X , ' CONST  t  WH1,DSLW,NT,R,DIBK,SH1,AIR  ?*,/) 

102  FORMAT (/, 3X, 'CONTINUB?  ("YBS"  OR  •'NO*')*,/) 

103  FORMAT (Al) 

104  FORMAT C'1',//,5X,'J.A.  WBISS,  *, 

2  'ADVANCES  IN  MICROWAVBS,  VOL.  8,',/7X,'PP.  295-320 j  *, 

3  'ACADEMIC  PRESS,  1974. ' ,// ,10X, 'PARAMETERS  OP  MICROSTRIP  *, 

4  'TRANSMISSION  LINBS' ,/,12X, 'AND  OP  COUPLED  PAIRS  OP  ', 

5  'MICROSTRIP  LINES') 

105  PORMATC// ,11X, 'H2/H1  -  ' ,P6 . 3 , 5X, *K  -  ' ,P6 . 3 , 4X, 'SINGLE  STRIP') 

106  FORMAT ( / , 10X , 'W/Hl* ,12X, 'ZO' ,14X, 'V' ,12X, 'K-BPP* ,11X, 'C' ,/, 

2  10X , 15X , ' OHMS ' , 8X , ' B*08  M/SBC' ,23X, 'PP/M' , /) 

107  FORMAT (//3X,'H2/H1  »  ',P6.3,4X,'R  -  ' , F6 . 3 , 4X , ' COUPLED  STRIPS', 

2  5X, 'S/HI  -  ' ,P5.2) 

108  PORMAT(/ ,3X, 'W/Hl' , 6X, ' ZO (B) ' ,4X, 'ZO(O) ' ,6X, ' V(B) ' ,4X, 'V(O) ' ,4X, 

2  'K-BPP (B)  K-BPP (0)',4X,'C(B)',5X,'C(0)', /18X, 'OHMS' ,12X, 

3  'B*08  M/SBC', 30X, 'PP/M',/) 

109  PORMATC//,'  AIR  SHOULD  BB  0.0  OR  ♦1.0',//) 

110  FORMAT ( / / , 5X , 2 (A8 , 2X) , / ) 

199  PORMATC///) 

901  CONTINUB 
999  CONTINUE 
CALL  EXIT 
STOP 
END 
C 
C 

SUBROUTINE  OUTPUT (CAP1B ,  CAPKB , CAPIO , CAPKO ,WH , IW ,NN) 

C 

IMPLICIT  RRAL*R (A-H.O-Z ' 

DATA  C/2.99792458/ 


MST01570 

MST01580 

MST01590 

MST01600 

MST01610 

MST01620 

MST01630 

MST01640 

MST01650 

MST01660 

M6T01670 

MST01680 

M6T01690 

MST01700 

MST01710 

MST01720 

MST01730 

MST01740 

MSTOWS? 

MST01760 

MST01770 

M6T01780 

MST01790 

MST0180C 

MST01610 

M6T018.20 

MST01830 

M6T0184Q 

MST01850 

MST01867 

MST01870 

MST01880 

MST01890 

MST01900 

MST01910 

MST01920 

H6T01930 

MST01940 

MST01950 

M5T01960 

MST01970 

MST01980 

MST01990 

MST02000 

MST02010 
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BFFKE-CAPKB/CAP1B 
RKE—DSQRT ( BFFKB ) 

ZOB-1 . QE*04 / (C*CAP1B*RKB) 

VELE«(1.0/RKE)*C 

IF (NN.  SO.  1)  WRITS ( IW,  101)  WH.ZOE, VBLB, BFFKB, CAPKB 
IF (UN  .SO.  1)  GO  TO  3 
EFFKO-CAPKO / CAPIO 
RKO-DSQRT(BFFKO) 

ZOO-1. 0B+04/ (C*CAP10*RK0) 

VBLO— ( 1 . D / RKO ) *C 

WRITB ( IW, 103 )  WH , ZOB , ZOO , VBL3 , VBLO , BFFKB , BFFKO , CAPKB , CAPKO 

101  FORMAT (5 ( 5X, FIG . 3) ) 

102  FORMAT (2X, F6 . 3 , 4X, F7 . 3 , 1X,P7. 3 , 4X,P6 . 3 , 2X, F6 . 3 , 5X, F6 . 3 , 2X, F6 . 3) 

103  FORMAT (2X,F6 . 3 , 4X,F7. 3 , IX, P7 , 3 , 4X, F6 . J  ,2.\,F6 .  3 , 5X,  P6 . 3 , 2X ,  F6 . 3 , 

2  4X,F8. 3 ,1X,F8 . 3) 

3  RBTORN 
BND 


S0BROCTI1JK  XGKS(M,WH,SH1) 

IMPLICIT  RBAL«8(A-H,0-Z) 
CO!*ION/DIM/X(5:>') , PHI  (59)  ,A(210)  ,B(210) 
AK-M 

WHM-WH/AM 

801  DO  901  I-l.M 
AI-I 

901  X(I)«(AI-1.0)»WHM 

IF (SHI  .80.  0.0)  GO  TO  1 
2,M-1 

802  DO  902  I-1,M< 

AI-I 

902  X  ( I +M) —At ‘WHM+SH1 


RBTORN 

KID 


SOBROOTINB  AMATfS.M) 

IHFLICIT  REAL*8 (A-H,0-Z) 
GCWMON/DIM/X(59) ,PHI(59) ,A(210) ,B(210) 

803  DO  903  I— 1,M 

INDBX1-M»1-I 

804  DO  904  J— I , INDBX1 
NF1-0 

805  DO  905  K«1,J 

905  NF1-HF1+K 
NF2-0 

IF (I  .LB.  2)  GO  TO  10 
INDBX2-I-2 

806  DO  906  L— 1 , INDHX2 

906  NF2— NF2*L 
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MST02020 

MST02030 

MST02040 

MST02050 

MST02060 

MTT02070 

MST02080 

MST02090 

M6T02100 

MST02110 

MST02120 

MST02130 

MST02140 

MST02150 

MST02160 

MCT02170 

M8T02180 

MST02190 

M8T02200 

MST02210 

MET02220 

MST02230 

MST02240 

M6T02250 

M6T02260 

M6T02270 

M6T02280 

M6T02290 

M6T02300 

M6T02310 

MST02322 

M6T02330 

MGT02340 

NWT02350 

M6T02360 

M6T02370 

MST02380 

MST02390 

M6T02400 

M6T02410 

M6T02420 

MST02430 

MST02440 

MGT02450 

MST02460 

MST02470 

MST02480 

,JST02490 

MST02500 

MSTC2510 

MST025Z0 

MST02530 


}  V  V-J  X  v 


10  N”NF1*NF2+(I-1)*J 

INDEX3-3*M»2*  (l-J)-I 
A(N)=PHI(I)*S*PHI(INDSX3) 

B(N)=PHI(I)--S*PHI(INDBX3) 

904  CONTINUE 
903  CONTINUB 
RETURN 
BHD 
C 
C 

SUBROUTINE  MPHI (NH,M,S) 

C 

IMPLICIT  REAL* 8 (A-H.O-Z) 

COWON/DIM/X(59j  ,PHI(59) ,A(210) ,B{210) 

AM-M 

INDBX1«3*M-1 

TP(S  .BO.  0.0)  IHDBX1-M 

8XW0N-WE/(2.0*AM) 

WXWON-l.OD  00 

807  DC  907  K-l.INDBXl 
XO-X(K) 

EXP— XO  ♦ EXWON 
BXN-XO-KXWON 
WYP— 2 . 0  *WYWON 

PHI (K)« (BXN/ (2. 0*KXWON) ) *DLOG( (EXN»*2) / (BXN*»2-M«P**2) ) - (BXP 

2  /  (2.0*KXW0N)  )  *DLOG(  (BXP**2) / (8XP*»2*WYP»»2)  )  ♦  (KYP/BXWON) 

3  *  (DATAN  (  HXP/HYP )  -DATAN  ( SXN/NXT  )  ) 

907  CONTINUB 

RETURN 

BND 

C 

C 

SUBROUTINE  M3RBBN{M,NH,SH1,DISK(R) 

C 

IMPLICIT  REAL'S (A-H.O-Z) 

COWON/DIM/X(59) ,PHI(59) ,A(210) , B(210) 

COWON  /PASS/  Rl,COfBO,DISKl 

EXTERNAL  GIST 

AM-M 

Rl-R 

DISX1— DISK 
CO-KH/AM'0 . 5 
XI— 5 . 0 
INT— 2 

H— Xl/DPLOAT ( 1ST) 

IP (SHI  .BQ.  0.0)  MA*M 
IP (SHI  ,NB.  0.0)  MA-3*M-1 

808  DO  908  t+t-l,MA 
BO— X  (W) 

YTOT-O . ODO 
XU-0.0 

XL— 0 , 0 


MST02540 

MST02550 

MST02560 

MST02570 

MST02580 

MST02590 

MST02600 

MST02610 

MST02620 

MST02630 

MST02640 

MST02650 

MST02660 

MST02670 

MST02680 

MST02690 

MST02700 

MST02710 

MST02720 

M6T02730 

MST02740 

MST02750 

MST02760 

MST02770 

MST02780 

MST02790 

MST02800 

MST02810 

M3T02320 

MSTQ2830 

MST02840 

MGT02850 

M6T02860 

M6T02870 

MST02880 

MST02890 

MST02900 

MST02910 

MST02920 

MST02930 

MST02940 

MST02950 

MST02960 

MST02970 

MST02980 

MST02990 

MST03000 

MST03010 

MST03020 

MST03030 

MST03040 

MST03050 


C 

C  COMPOTE  FIRST  INTEGRAL 

C 

809  DO  909  I»1,INT,1 
XU=*XU*H 

CALL  DGG32(XL,XU,GIHT,YRSSUL) 
VTOT«VfT(jT*  YRBSUL 


MST03060 

MST0307G 

MST03080 

MST03090 

MST03100 

MST03110 

KST03120 


XL»XL*H 
909  CONTINUE 
AI1«»YTGT 
C 

C  COMPUTB  SECOND  INTEGRAL 

C 

Sl»(CO+BO)«Xl 

32«<B0-CC)*X1 

CALL  SICI(SI,CI.-S1) 

AI2A«DSXN<Sl)/m-(CO*BO;*CI 

CALL  SICI(SI,CIrS2) 

AI2B~DSIH (S2) /XI- (BO-CO) *CI 
AI2-AI7A-AI2B 

PHI  (JK)»>4.G*  (All *1.0/  ( (1. 0+DIBK)  *CO*2.0)  *AI2) 
908  CONTINUE 
RETURN 
END 

FUNCTION  GINT  (U) 

IMPLICIT  R3AL*8(A-!i,0-2) 

COWON  /PASS/  R, CO, BO, DIRK 

V-(R-1.0)*U 

Wl-CO»U 

H2-DC08H(U) 

W3-BSINH(U) 

W4-DCOS (BO*U) 

IF (R  .NS.  0)W4*W4*D6IirHC?) 

IF(R  .EC.  0)DEN«W3<-DIER*W2 

IF (R  .HH.  0)DHN"W3*DCOSH (V) ♦DISK*W2*DSINK(V) 

GINT-DS  IN  ( W1 )  /HI  *W3  /  C  »N4  /  DEN 

RETURN 

BHD 


KST03130 

HST03140 

M5TG3130 

KST03160 

MST03170 

MST03130 

MST03190 

HST03200 

MST03210 

MST03220 

MST03230 

MST03240 

MST03250 

Mb‘5’33260 

MST03270 

MET03230 

MST03290 

MST03300 

HST03310 

M6'r0332Q 

MST03330 

MST03340 

MET03350 

MST0334D 

H6T03370 

MST03380 

KCT03390 

MST03400 

MST03410 

MST03420 

MST03430 
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